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ABSTRACT 


The potential for o:ean prediction within the next decade is discussed in terms 
of the crucial elements: data, computing power, and models. The paper outlines a 
strategy for a global ocean prediction system that is based on appropriate satellite 
data, Class 7 computers (~1 gigaflop and 32-128 million words), and eddy-resolving 
basin or global scale models. The paper is focused on the prediction of meandering 
currents, frontal positions, and eddies, but prediction of other oceanic phenomena 
is also surveyed. 

1. INTRODUCTION 

The next decade should be to numerical prediction of the ocean circulation 
what the 1950’s and 1960’s were to numerical weather prediction, a time when the 
essential elements came together to permit major advances in forecasting skill. 

These elements are 1) adequate data input, 2) adequate computing power, and 3) 
properly designed and adequately validated models for both data assimilation and 
forecasting. Ten vears ago none of these requirements were satisfied for many 
aspects of ocean prediction. Meandering currents and eddies provide a notable 
example: (1) clear prospects for adequate data with high density and global coverage 
from satellite-borne instruments have become evident only in the last few years, 

(2) Class 7 computers with a sustainable speed of ""1 gigaflop and 32-128 M words of 
memory are required for eddv-resolving forecasts in major ocean basins, and (3) 
notable differences from meteorology are required in model design and in data 
assimilation. 

The U.S. Navy has determined a requirement for an ocean prediction capability, 
and has established a research program to meet the requirement. This paper is an 
overview of the requirements for data and computing power, and the issues of model 
design and data assimilation primarily in the context of the effort at the Naval 
Ocean Research and Development Activity. It should be stressed that there are many 
types of ocean prediction aimed at different aspects of the motion, internal 
structure, and surface of the ocean. These provide several important roles for 
satellite data. This paper is focused on one important aspect of ocean prediction 
where altimeter data can play a central role. In particular, there will be a 
demonstration of the role that satellite altimeter data can play in detailed 
surface and subsurface forecasts of meandering currents, eddies, and frontal 
locations, forecasts that may extend up to several months. 

2. DIFFERENT TYPES OF OCEAN PREDICTION 

2 . 1 Classes of Response to Atmospheric Forcing 

Before concentrating on a particular type of ocean prediction, it is useful 
to survey the scope of the problem. Table 1 from Hurlburt (1984) shows three 
clas ses of oceanic response to atmospheric forcing that cover many phenomena for 
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♦From Hurlburt (1984) 



which prediction appears to be feasible and useful. It is not all-encompass ing 
(e.g. tides and tsunamis are omitted), but it is Intended to help put the multi- 
faceted problem of ocean prediction into perspective, identify phenomena where 
predictive skill is feasible, clarify the potential roles of satellite data, and 
to define the focus of the paper, which is primarily on Class 2. Fortunatelv, 
many phenomena are sufficiently decoupled that it is possible to design forecast 
models that can simulate and predict certain phenomena while suppressing or 
parameterizing the statistical effects of the remainder. 

2 . 2 No wcas ting 

"Nowcasting 1 ’ is important in all three classes. To nowcast the models are 
integrated forward in time while driven by new atmospheric forcing functions and 
assimilating appropriate new ocean data as they become available. This allows 
the models t< 

a. Fill in temporal gaps in the data by using theLr predictive skill. 

b* Convert better observed surface fields into subsurface structure, 

c. Convert better observed atmospheric forcing functions into oceanic 
information . 

2 . 3 Class 1: Strong, Rapid, and Direct 

Until now ocean prediction efforts have concentrated primarily on tides and 
certain phenomena in Class 1, a class which to some extent ir^ludes the cumulative 
effects of fine scale phenomena that are predicted via empirical or semi-empxrical 
parameterizations. Storm surges have received the moat attention (Welander, 1961; 
Jelesnianski, 1967; Crawford, 1979), but recently the U.S. Navy has initiated 
hemispheric forecasts of surface waves (Pierson, i982) and the upper mixed layer 
(Clancy and Martin, 1981; Clancy and Poliak, 1°83). 

The output from atmospheric ir, and satellite-borne instruments 

measuring surface wind speed and dirowdon (^c^cercmeter) or even wind speed alone 
(scanning multi-frequency microwave rad (SMMR)) show great promise in 
facilitating predictions in Class 1 and cl^ss 3. Scatterometers and microwave 
radiometers can provide adequate rsvera^e, accuracy, and spatial and temporal reso- 
lution except for coverage boundaries (Jones et al, 1982; Wentz et al, 

1982; Lipes, 1982; Satellite Surface Stress Working Group, 1982; 'Mueller, 1982). 

For mixed layer forecasting, measurements of sea surface temperature (SST) (multi- 
channel infrared and microwave radiometers) would be of value. Some work has also 
been done or satellite measurement of latent heat flux (Liu, 1984) and Incident 
solar radiation at the surface (Gautier, 1981), the two heat fluxes that are usually 
the largest at the air-sea interface. 

Although a satellite altimeter can measure surface wind speed and significant 
wave height in addition to sea surface elevation (Fedor and Brown, 1982), the value 
of a single altimeter with only a nadir beam is severely degraded when applied to 
prediction of phenomena in Class 1. This is due to the inverse relation between 
spatial and temporal resolution that prevents adequate resolution for these pheno- 
mena in both space and time simultaneously. However, the altimeter can be useful 
in constructing climatologies of wind speed and significant wave height (Chelton 
et al., 1981; Wentz et al., 1982; Mognard et al., 1983), in validation and tuning of 
certain models and in some research problems involving these quantities. Because 
tides are repetitive with known periods, the altimeter should be important in 
refining the global knowledge of tidal phase and amplitude (Cartwright and Alcock, 
1981; Brown and Hutchinson, 1981; Diamante and Nee, 1981; Mazzega, 1983). 
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2.4 Class 2: Slow and Indirect 


The prediction of mesoscale (~ 50-500 km) eddies, meandering currents, and 
frontal positions (phenomena in Class 2) is the focus of this paper. It poses a 
problem where altimeter data should play a central role. This Is because (1) 
individual mesoscale eddies and current meanders are often not driven direct y by 
the wind or by any other external forcing function, and (?) eddies can persist for 
more than a year after their initial generation (Lai and Richardson, 1977). Thus, 
oceanic data are crucial for reliable prediction of individual eddies and current 
meanders. With sufficient oceanic data, we anticipate that the prediction of these 
features can be treated as an initial value problem in which the future forcing 
functions are representative but not accurate on time scales greater than & few days. 
Without oceanic data input, simulation can he used to predict the statistical prop- 
erties of features and ensembles, but usually not the evolution and movement of 
individual features. 

There has been little work on the prediction of Class 2 phenomena except for 
the work on prediction of mesoscale eddies, using limited-area models, by A.R. 
Robinson's group at Harvard (Robinson and Haiavogel, 1980; Robinson and Tu, 1982; 
Miller and Robinson, 1983). However, there is a substantial body of literature on 
the simulation and ocean dynamics of Class 2 phenomena which is very useful in the 
design of ocean circulation prediction models (Holland and Lin, 1975; Rhines, 1977; 
Semtner and Mintz, 1977; Holland, 1978, 1982; McWilliams and Flierl, 1979; Cox, 

1979; Hurlburt and Thompson, 1980, 1982; Lin and Hurlburt, 1981; Schmitz and 
Holland, 1982; Heburn et al. , 1982). Many of these have simulated some features of 
the ocean circulation with substantial success using simple models, rectangular 
domains, and simple steady forcing functions. This is an indication of potential 
predictive skill, but one largely unverfied at this time. 

As an example. Figure 1 illustrates model simulation skill for a Class 2 
phenomenon in the Gulf of Mexico. The model has two active layers, realistic 
geometry and bottom topography, and horizontal resolution sufficient for major 
current systems and eddies. Gurnrisingly , no eddy-resolving ocean modeling study 
published to date has included all of these features in a deep ocean basin where 
planetary vorticity advection va.> important and the model was integrated to 
statistical equilibrium. This calculation was performed by A. Wallcraft (personal 
communication) using the model of Hurlburt and Thompson (1980). The ability to 
handle detailed coastline geometry was added by Wallcraft. The model was driven 
from rest to statistical equilibrium solely by a steady inflow through the Yucatan 
Straits between Cuba and Mexico which was compensated by outflow through the 
Florida Straits. 

Figure 1 compares "instantaneous" upper ocean flow patterns (a) from the 
numerical model and (b) from observations by Leipper (1970). The Loop Currant is 
the major current system depicted. At this point it is about to shed an eddy. 

The Loop Current is observed to penetrate into the Gulf, bend westward, and to 
shed an eddy with a period of approximately one year. T is "annual" cycle was 
„ long thought to be due to seasonal variations in the flow through the Yucatan 

Straits (Cochrane, 1965). However, The model Loop Current exhibits an approximately 
annual eddy-shedding period when the Inflow is steady, contradicting the earlier 
hypothesis. Although time variations are not essential, they can still play a 
significant role in the eddy shedding (Hurlburt and Thompson, 1980). The model 
7" » Loop Current also spontaneously shed eddies with realistic diameters, amplitudes, 
and westward propagation speeds. The dynamical basis for the agreement between the 
IV observations and a circulation model with low vertical resolution la discussed by 
Hurlburt and Thompson (1980, 1982). Additional information on the simulation in 
Figure 1 is found in Hurlburt (1984), including the spontaneous generation of a 
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Figure 1. (a) Instantaneous view of the interface deviation in a two-layer 

simulation of the Gulf of Mexico driven from rest to statistical equilibrium 
solely by inflow through the Yucatan Straits. The contour interval is 25 m 
with solid contours representing downward deviations, (b) Depth of the 22°C 
isothermal surface, 4-18 Aug 1966 (Aliminos cruise 66-A-ll) from Leipper 
(1970). The contour interval is 25 m. (From Hurlburt, 1984) 
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counter-rotating vortex pair as i;he eddy propagates westward. Such a feature has 
been observed repeatedly in the western Gulf of Mexico. 

2.5 Class 3: Slow and Direct 


Phenomena in Class 3 are a direct response to atmospheric forcing, but on much 
longer time scales than those in Class 1. Because they are a much more integrated 
response to the wind, the requirements for temporal resolution and accurate depic- 
tion of daily atmospheric fluctuations are not as stringent as in Class 1. This 
should permit forecasts on time scales greater than the few-day time scale for 
atmospheric predictability. In Class 3 "nowcasting" and forecasting may be 
feasible by using new atmospheric forcing functions without new oceanic data. 

The addition of altimeter-derived sea surface elevation accurate on the scales 
of maior ocean basins as well as on mesoscales would reduce the burden on the 
simulation skill of the model. Examples of Class 3 simulation and prediction 
studies include the work by Haney (1980) on predicting large-scale anomalies in 
the North Pacific, by Busalacchi et al (1983) on the tropical Pacific, and by 
Preller and Heburn (1983) on the western Mediterranean Sea. 

3. RESOLUTION REQUIREMENTS FOR ALTIMETER AND SCATTEROMETER DATA 

The horizontal and temporal resolution required for the forcing functions is 
dictated primarily by the space and time scales of the atmosphere. Thus, scatter- 
ometer -derived wind stress data would be needed on a daily basis to resolve the 
evolution of atmospheric storm systems. The spatial resolution of the wind stress 
data should be sufficient ("' 50 km) to resolve the wind stress curl associated with 
atmospheric cyclones, anti-cyclones and fronts. V'nds as close to coastlines as 
possible are highly desirable. In these regions horizontal resolution should be 
2C km or better. Since a satellite scatterometer with an orbit altitude of about 
850 km has about a 1500 km swath width (with gaps) , a single satellite in conjunction 
with atmospheric models and the existing data base could meet these requirements 
adequately (Satellite Surface Stress Working Group, 1982). 

Altimeter data is essential for global prediction and representation of indivi- 
dual current meanders and eddies. Cheney and Marsh (1981) have already demonstrated 
the ability of the SEASAT altimeter to detect the Gulf Stream and mesoscale eddies 
when using either repeat tracks or a geoid, and Cheney et al (1983) have used nine 
sets of 3-day repeat tracks from SEASAT to produce a global map depicting the RMS 
mesoscale variability of the sea surface elevation during Sep - Oct, 1978. Tapley 
et al (1982) and Wunsch and Gaposhkin (1980) review the problem of obtaining 
accurate measurements of current-related variations in the sea surface elevation 
from a satellite altimeter and the success achieved with SEASAT. To be useful in 
ocean circulation prediction, the altimeter data must, of course, resolve major 
current systems and eddies in time, amplitude, and horizontal dimension. However, 
the amplitude dependence on horizontal scale should also be noted. For currents 
in geostrophic balance 

k x f Vg - -gVn 

where f is the Coriolis parameter, Vg is the velocity of the geostrophic surface 
current, g is the acceleration due to gravity, and n is the sea surface elevation 
above the geoid related to oceanic surface currents. The maximum velocity (V) 
for typical, significant, persistent, mesoscale (~ 50 to 500 km) features in the 
ocean ranges from about 10 cm/sec to 1 m/sec. Up to 2 m/sec is not uncommon in 
the stronger currents. For an altimeter which can usefully measure oceanographic 
features with a minimum amplitude of r) * 10 cm, the minimum radius of an eddy or 
half-width of a current (L) with V ■ 10 cm/sec is L « 100 km at 42° Lat and L ■ 150 
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km at 27° Lat. For V = 100 cm/sec, L - 10 km at 42° Lat , and the Rossby number, 
Rq = V/fL = 1. For R 0 « 1 the flow is geostrophically balanced. For R c >> 1 
vortices can develop cyclostrophic balance. 


2 



where vq is the tangential velocity and r is the radius of curvature. In this 
case the amplitude of n associated with the vortex is dependent on v^ but is 
nearly independent of the scale of the vortex. However, a velocity scale > 2 m/s 
is required for a signal in the sea surface elevation > 10 cm. 

From this we conclude that a satellite altimeter able to measure variations 
in the sea surface elevation > 5 cm, which are associated with oceanic currents, 
should provide usable information for the more significant currents and eddies 
with scales greater than a few tens of kilometers. This accuracy is feasible 
with current technology (TOPEX Science Working Group, 1981). Altimeter track 
spacing of about 30 km and at most 50 km at mid-latitudes is needed for horizontal 
resolution of most detectable eddies. Experiments by J. Kindle at NORDA have 
shown that three ascending or descending tracks with error-free data are sufficient 
to adequately map an eddy with a nearly circular shape. Repeat intervals up to a 
month should provide adequate temporal resolution for most mesoscale features, 
particularly with the aid of predictive models to fill in the temporal gaps. 

To minimize the time required for altimeter tracks to span a given eddy, the orbit 
should be chosen so that adjacent tracks are one day apart. This requirement 
should override any concern about tidal aliasing, since the tides are relatively 
large scale and because removal of the tidal contributions to the sea surface 
elevation appears feasible (Schwiderski, 1980; Cartwright and Alcock, 1981; Brown 
and Hutchinson, 1981; Diamante and Nee, 1981; Mazzega, 1983). However, as noted 
in some of the preceding references, accurate determination of the tidal phases 
and amplitudes would be facilitated by at least once having an altimeter in an 
orbit relatively free of tidal aliases with long periods. 

A single satellite in a 20-day repeat orbit which carried a multi-beam alti- 
meter (Bush et al, 1984) could easily meet the preceding requirements. Minimal 
satisfaction of the requirements could be achieved by using a single satellite 
with a nadir beam altimeter in a 40-day repeat orbit, if a numerical model was 
able to bridge the temporal gaps. In the absence of an adequate geoid, it would 
be necessary to use a mean sea surface, so that only deviations from the mean would 
be available to the model. In that case the model could obtain the mean from its 
own or observed climatology. Bandpass filtering could be required to remove short 
(< 0(10 km)) and long wave (> a few thousand km) errors, the scales with the 
greatest errors (TOPEX Science Working Group, 1981). 

4. STRATEGY FOR OCEAN PREDICTION 

Highly efficient models with low vertical resolution have demonstrated a 
remarkable ability to simulate meandering currents and eddies, but high horizontal 
resolution (■* 10 km) is required to resolve them. To avoid serious limitations on 
data assimilation and predictive skill, it is advantageous to use domains which 
cover major ocean basins or serai-enclosed seas where the phenomena of interest 
are generated within the domain, and where the open boundary segments are either 
small (such as flow through a strait) or very distant from the region of interest. 
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Circulation models on subdomains of maior ocean basins which have extensive open 
ocean boundaries are extremely dependent on -^e source like a large-scale model 
for boundary conditions during both the data assimilation and prediction phases 
(Miller and Robinson, 1983). To provide useful boundary conditions, the large- 
scale model must also have resolution sufficient for the currents and eddies. To 
model major ocean basins with high horizontal resolution and low vertical resolution, 
a Class 7 computer 1 Gigaflop and 32M to 128M words) is required. The intro- 
duction of Class 7 computers is anticipated in the middle and late 1930* s. 

Cur^ntly, satellites and atmospheric forecast models provide the only real 
prospects for oceanic data and forcing functions with global coverage and resolution 
adequate for use in ocean circulation models. For this purpose, the fields of 
greatest value, which are observable from satellites, are altimeter-derived sea 
surface elevation, scatterometer-derived wind stress, and heat fluxes. Successful 
ocean prediction appears feasible without extensive subsurface data collection by 
using the simulation skill of the model to convert the potentially well-observed 
fields at the surface into subsurface information. This looks more promising 
than attempting to assimilate the extremely sparse subsurface data, a process 
likely to do more harm than good. Where eddy-resolving subsurface data with 
useful coverage is available (nowhere at present), it should, of course, be used, 
especially if it is available on a regular basis and its statistical properties 
are well known. In the long term, acoustic tomography shows promise in this 
area (Munk and Wunsch, 1982), but we anticinate that subsurface data acnuisition 
adequate for prediction of mesoscale features will, at best, be limite in 
coverage for at least the next decade. Some subsurface data and major field 
programs are essential for the local and regional subsurface validation of the 
forecast models. 

The surface to subsurface data conversion process would be accomplished bv 
using both the altimeter data and the atmospheric data as forcing functions. 

Although this is an untested hypothesis at this time and the assimilation of the 
altimeter data could be troublesome, the potential for surface to subsurface 
conversion is enhanced by noting that (1) a single internal vertical mode explains 
much of the oceanic variability (Richman et al, 1977; Flierl, 1978), (2) in a 
layered circulation model with a single internal mode there is a one-to-one 
correspondence between variations in the sea surface elevation and variations in 
the depth of the pycnocline, (3) ©uch a model can be initialized by altimeter 
data alone, except near the equator where wind stress is necessary, and (4) such 
models have demonstrated remarkable simulation skill in certain situations 
(Busalacchi et al, 1983; Lin and Hurlburt, 1981; Hurlburt and Thompson, 1980, 1982), 
although they miss the sometimes important effects of topography and baroclinic 
instability (Holland and Lin, 1975; McWilliams et al, 1978; Hurlburt and Thompson, 
1983). These require at least one additional vertical mode and we recommend no 
less than three in an operational forecast model. High vertical resolution 
can be obtained by coupling the circulation model with low vertical resolution to 
an array of one-dimensional mixed layer models with high vertical resolution. A 
hemispheric array makes up the mixed-layer forecasts by the U. S. Navy (Clancy 
and Martin, 1981). 

Among the ideas proposed for application to ocean prediction, some that show 
promise are 1) the use of satellite altimeter and scatterometer data, 2) the 
conversion of surface into subsurface information, 3) the partial decoupling 
of the ocean circulation problem from the turbulent thermodynamic mixed layer 
problem, 4) eddy-resolving circulation models with low vertical resolution and 
basin or global-scale domains, and 5) Class 7 computers. Some that are not 
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promising (at least within the next decade or more) are l) eddy-resol ving models 
for major ocean basins with many fixed levels in the vertical, 2) the inclusion 
of sophisticated mixed layer physics in eddy-resolving circulation models 
covering major ccean basins, 3) the assimilation of sparse in situ data by ocean 
circulation models, 4) coarse grid, non-eddy-resolving prediction models for the 
ocean circulation, and 5) sub-basin-sca! models with extensive open boundaries 
except when they can obtain boundary conditions from an eddy-resolving basin 
scale model. "Nowcasts" and short-range forecasts (a few days) from stand-alone 
limited-area models appear feasible in regions with adequate subsurface sampling 
(none at present). The first two ideas lack promise only because of the anti- 
cipated computing power within the next decade or more. 

Within the next decade the appropriate satellite data, Class 7 computers, and 
eddy-resolving basin or global-scale models should become available, and they should 
be used to form the heart of a global ocean prediction system. In the meantime, 
alternate methods can provide a more limited capability which is applicable in 
some situations. Forcing functions from atmospheric models should permit useful 
nowcasting and predictive skill for Class 1 and Class 3 phenomena listed in Table 
1 prior to the availability of satellite altimeter and scatterometer data. 

ACKNOWLEDGEMENTS 

Funding is from the Office of Naval Research under the Special Focus Program, 
"Ocean dynamics from OEOSAT". This paper is a condensation of Hurlburt (1984). 

The manuscript was typed by Charlene Parker. 




REFERENCES 




Brown, R. D,, and M. K. Hutchinson, 1981: Ocean tide determination from satellite 

altimetry. Oceanography from Space . J. F. R. Gower, Ed. Plenum Press, 

New York, 897-906. 

Busalacchi, A. J., K. Takeuchi, and J. J. 0 ! Brien, 1983: On the interannual 

wind-driven response of the tropical Pacific Ocean. Hydrodynamics of the 
Equatorial Ocean , J. C. J. Nihoul, Ed., Elsevier Science Publishers, 
Amsterdam, 155-195. 

Bush, G. B., E. B. Dobson, R. Matyskiela, C. C. Kilgus, and E. J. Walsh, 1984: 

An analysis of a multibeam altimeter. Marine Geodesy (in press). 

Cartwright, D. E. and G. A. Alcock, 1981: On the precision of sea surface 

elevation and slopes from SEASAT altimetry of the northeast Atlantic Ocean. 
Oceanography from Space , J. F. R. Gower, Ed., Plenum Press, New York, 

885-895. 


Chelton, D. B., K. J. Hussey, and M. E. Parke, 1981: Global satellite measure- 

ments of water vapour, wind speed and wave height. Nature , 294 , 529-532. 

Cheney, R. E. and J. G. Marsh, 1981: SEASAT altimeter observations of dynamic 

topography in the Gulf Stream region. J. Geophys. Res. , 86 , 473-483. 

Cheney, R. E., J. G. Marsh, and B. D. Beckley, 1983: Global mesoscale variability 

from collinear tracks of SEASAT altimeter data. J. Geophys. Res., 88, 
4343-4354. 


Clancy, R. M. and P. J. Martin, 1981: Synoptic forecasting of the oceanic 

mixed layer using the Navy’s operational environmental data base: 
present capabilities and future applications. Bull. Am. Meteorol. Soc., 
62, 770-784. 

Clancy, R. M. and K. D. Poliak, 1983: A real-time synoptic ocean thermal 

analysis/forecast system. Prog. Oceanog. , 12 , 383-424. 

Cochrane, J. D., 1965: The Yucatan Current and equatorial currents of the 

western Atlantic. Unpublished report, Dept, of Oceanography, Texas A&M 
University, Ref. (65-17T), 20-27. 


Cox, M. D., 1979: A numerical study of Somali Current eddies. J. Phys. Oceanogr 

9, 311-326. 


Crawford, K. C., 1979: Hurricane surge potentials over southeast Louisiana as 

revealed by a storm-surge forecast model: a preliminary study. Bull. Am. 

Met. Soc. , -60, 422-429. 



Diamante, J. M. and T. S. Nee, 1981: Application of satellite radar altimeter 

data to the determination of regional tidal constituents and the mean sea 
surface. Oceanography from Space , J. F, R. Gower, Ed., Plenum Press, 

New York, 907-918. 

Fedor, L. S. and G. S. Brown, 1982: Waveheight and wiri speed measurements from 

the SEASAT radar altimeter. J. Geophys. Res. , 87 , 3254-3260. 


596 




Flierl, G. R., 1978: Models of vertical structure and the calibration of 

two-layer models. Dyn. Atmos. Oceans , 2^ 341-381. 

Gautier, C., 1981: Daily shortwave energy budget over the ocean from geo- 

stationary satellite measurements. Oceanography from Space , J. F. R. Gower, 
Ed., Plenum Press, New York, 201-206. 

Haney, R. L. , 1980: A numerical case study of the development of large-scale 

thermal anomalies in the Central North Pacific Ocean. J. Phys. Oceanogr. , 
10 , 541-356. 

Heburn, G. W., T. H. Kinder, J. H. Allender, and H. E. Hurlburt, 1982: A 

numerical model of eddy generation in the southeastern Caribbean Sea. 
Hydrodynamics of Semi-enclosed Seas , J. C. J. Nihoul, Ed., Elsevier 
Scientific Publishing Company, Amsterdam, 299-327. 

Holland. W. R., 1978: The role of mesoscale eddies in the general circulation 

of the ocean - numerical experiments using a wind-driven quasi-geostrophic 
model. J . Phys . Oceanogr . , _8, 363-392. 

Holland, W. R., 1982: Simulation of North Atlantic variability. Proceedings 

of the Workshop on Gulf Stream Structure and Variability , J. M. Bane, Jr., 
Ed., Sponsored by U. S. Office of Naval Research, 114-126. 

Holland, W. R. and L. B. Lin, 1975: On the generation of mesoscale eddies and 
their contribution to the oceanic general circulation. J . Phys . Oceanogr . , 
5, 642-669. 

Hurlburt, H. E., 1984: The potential for ocean prediction and the role of 

altimeter data. Marine Geodesy (in press). 

Hurlburt, H. E. and J. D. Thompson, 1980: A numerical study of Loop Current 

intrusions and eddy shedding. J . Phys . Oceanogr . , 10 , 1611-1651. 

Hurlburt, H. E. and J. D. Thompson, 1982: The dynamics of the Loop Current 

and shed eddies in a numerical model of the Gulf of Mexico. Hydrodynamics 
of Semi-enclosed Se as, J. C. J. Nihoul, ed., Elsevier Scientific Publishing 
Company, Amsterdam, 243-297. 

Hurlburt, H. E. and J. D. Thompson, 1983: Preliminary results from a numerical 

study of the New England Seamount Chain influence on the Gulf Stream. 
Predictability of Fluid Motions , G. Holloway and B, J. West, Eds., American 
Institute of Physics, New York, 489-504. 

Jelesnianski, C. P., 1967: Numerical computations of storm surges with bottom 

stress. Mon . Wea . Rev . , 95 , 740-756. 

Jones, W. L., L. C. Schroeder, D, H. Boggs, E. M. Bracalente, R. A. Brown, 

G. J. Dame, W. J. Pierson, and F. J. Wentz, 1982: The SEASAT-A satellite 

scatterometer : the geophysical evaluation of remotely sensed wind vectors 

over the ocean. J. Geophys. Res. , 87 , 3297-3317. 

Lai, D. Y. and P. L. Richardson, 1977: Distribution and movement of Gulf Stream 

rings. J. Phys. Oceanogr. , 2* 670-683. 


597 



Leipper, D. F., 1970: A sequence of current patterns in the Gulf of Mexico. 

J. Geophys. Res. , 75 , 637-657. 

Lin, L. B. and H. E. Hurlburt, 1981: Maximum simplification of nonlinear 

Somali Current dynamics. Monsoon Dynamics , M. J. Lighthill and R. P. 

Pearce, eds., Cambridge University Press, 541-555. 

Lipes, R. G., 1982: Description of SEASAT radiometer status and results. 

J . Geophys . Res . , 87 , 3385-3395. 

Liu, W. T. , 1984: Estimation of latent heat flux with SEASAT SMMR, a case study 

in the North Atlantic. Large Scale Experiments and Satellites , C. Gautier 
and M. Fieux, Eds., Reidel Press, Boston (in press). 

Mazzega, P., 1983: The M2 oceanic tide recovered from SEASAT altimetry in 

the Indian Ocean. Nature , 302 , 514-516 

McWilliams, J. C. and G. R. Flierl, 1979: On the evolution of isolated, 

nonlinear vortices. J . Phys . Oceanogr . , 9^, 1155-1182. 

McWilliams, J. C., W. R. Holland, and J. H. S. Chow, 1978: A description of 

numerical Antarctic Circumpolar Currents. Dyn. Atmos. Oceans , 2 , 213-291. 

Miller, R. N. and A. R. Robinson, 1983: Dynamical forcast experiments with 

a baroclinic quasi geostrophic open ocean model. Predictability of Fluid 
Motions , G. Holloway and B. J. West, Eds., American Institute of Physics, 

New York, 457-472. 

Mognard, N. M., W. J. Campbell, R. E. Cheney, and J. G. Marsh, 1983: Southern 

ocean mean monthly waves and surface winds for winter 1978 by SEASAT radar 
altimeter. J. Geophys. Res., 88, 1736-1744. 

Mueller, J. L., 1982: Analysis of satellite remote sensing systems. Proceedings 

of the Ocean Prediction Workshop . C. N. K. Mooers, S. A. Piacsek, and A. R. 
Robinson, Eds. Sponsored by the Oceanographer of the Navy and the Office of 
Naval Research, 292-307. 

Munk, W. and C. Wunsch, 1982 : Observing the ocean in the 1990*8. Phil. Trans. 

R. Soc. Lond. . A, 307, 439-464. 

Pierson, W. J. 1982: Spectral Ocean Wave Model (SCWM) , a northern hemisphere 

computer model for specifying and forecasting ocean wave spectra. Tech. 
Report, DTNSRDC-82/011. David Taylor Naval Ship Research and Development 
Center, Bethesda, MD 186 pp. 

Preller, R. H. and G. W. Heburn, 1983: A simple numerical model of circulation 

in the western Mediterranean Sea. (Submitted to Oceanologlca Acta ) . 

Rhines, P. B., 1977: The dynamics of unsteady currents. The Sea , Vol. 6, 

E. D. Goldberg, I. N. Me Cave, J. J. O'Brien and J. H. Steele, Eds., Wiley 
Interscience, 189-318. 

Rlchman, J. G., C. Wunsch, and N. G. Hogg, 1977: Space and time scales of 

mesoscale motion in the western North Atlantic. Rev, of Geophys. and Space 
Phys. . 15 , 385 - 420 . 


598 






Robinson, A, R. and D, B. Haidvogel, 1980: Dynamical forecast experiments with 

a barotropic open ocean model. J . Phys . Oceanogr . , 10 , 1909-1928. 

Robinson, A. R. and K. Tu, 1982: A combined statistical/dynamical approach 

to regional forecast modeling of open ocean currents. Proceedings of the 
Ocean Prediction Workshop , C. N. K. Mooers, S. A. Piacsek, and A. R. Robinson, 
Eds. Sponsored by the Oceanographer of the Navy and the Office of Naval 
Research, 114-141. 

Satellite Surface Stress Working Group, 1982: Science Opportunities Using 

Satellite Wind Stress Measurements over the Ocean. J. J. O’Brien, Chairman. 
Prepared for: National Aeronautics and Space Administration, 153 pp. 

Schmitz, W. J., Jr. and W. R. Holland, 1982: A preliminary comparison of 

selected numerical eddy-resolving general circulation experiments with 
observations, J. Mar. Res. , 40 , 75-117. 

Schwiderski, E. W., 1980: On charting global tides. Rev, of Geophys. and 

Space Phys. , 18 , 243-268. 

Semtner, A. J., Jr. and Y. Mintz, 1977: Numerical simulation of the Gulf 

Stream and mid-ocean eddies, J. Phys. Oceanogr. , _7> 208-230. 

Tapley, B. D., G. H. Born, and M. E. Parke, 1982: The SEASAT altimeter data 

and its accuracy assessment. J. Geophys. Res. , 87 , 3179-3188. 

TOPEX Science Working Group, 1981: Satellite altimetric measurements of 

the ocean. C. Wunsch, Chairman. Report issued by NASA and JPL, 78 pp. 

Welander, P., 1961: Numerical prediction of storm surges, Advances in 

Geophysics , 8 , 315-379. 

Wentz, F. J., V. J. Cardone, and L. S. Fedor, 1982: Intercomparison of wind 

speeds inferred by the SASS, Altimeter and SMMR. J. Geophys. Res. , 87 , 
3378-3384. 

Wunsch, C. and E. M. Gaposhkin, 1980: On using satellite altimetry to determine 

the general circulation of the oceans with application to geo id improvement. 
Rev, of Geophys. and Space Phys. , 18, 725-745. 



